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ABSTRACT 

Recently, radial velocities have been measured for a large sample of M giants from the 2MASS catalog, selected 
to be part of the Sgr dwarf leading and trailing streams. Here we present a comparison of their kinematics to 
models of the Sgr dwarf debris orbiting Galactic potentials, with halo components of varying degrees of flattening 
and elongation. This comparison shows that the portion of the trailing stream mapped so far is dynamically young 
and hence does not provide very stringent constraints on the shape of the Galactic dark-matter halo. The leading 
stream, however, contains slightly older debris, and its kinematics provide for the first time direct evidence that 
the dark-matter halo of the Galaxy has a prolate shape with an average density axis ratio within the orbit of Sgr 
close to 5/3. 
Subject headings: dark-matter - Galaxy: halo, kinematics and dynamics, structure, fundamental parameters 



L INTRODUCTION 

Using M giant candidates from the 2MASS catalog, Majew- 
ski et al. (2003) have been able to trace streams from the Sagit- 
tarius dwarf through roughly 360 degrees on a great circle on 
the sky. More recently, the radial velocities of a subset of sev- 
eral hundred of these M giant candidates have been measured. 
The results have been presented in Law et al. (2003) and in 
tabular form (for the trailing stream) in Majewski et al. 2004 
(hereafter M04). Hence 4 of the 6 phase-space coordinates of 
the stars in the streams are now available for a comprehensive 
dynamical study. 

One of the intriguing questions that can in principle be ad- 
dressed with such a dataset, concerns the shape of the dark- 
matter halo of the Galaxy. Is this flattened? Prolate or oblate? 
Ibata et al. (2001) and Majewski et al. (2003) have argued 
that the fact that the debris of Sgr is distributed along an almost 
complete great circle on the sky implies that the dark halo of 
the Milky Way has to be very close to spherical, with allowed 
(density) axis ratios in the range 0.9 to 1 (Law et al. 2003). 
On the other hand, Helmi (2004) has shown that, if the debris 
discovered so far is made up of material lost in the last 3 to 4 
pericentric passages (i.e. less than 3 Gyr ago), then such debris 
would also be found along a great circle on the sky. Flattened 
(oblate) or elongated (prolate) halos with average density axis 
ratios, measured within the region probed by the debris, that are 
as low as 3:5 or as high as 5:3, respectively, would be allowed. 
Such values would be consistent with those found in numerical 
simulations of cold dark-matter halos (e.g. Bullock 2002). 

The motivation behind this Letter is to establish whether the 
additional kinematic information that has now become avail- 
able could perhaps break this degeneracy and allow us to better 
constrain the axis ratio of the Galaxy's dark-matter halo. 

2. COMPARISON OF MODELS AND DATA 

The models presented here are more thoroughly discussed 
in Helmi (2004). Now we only present a brief summary. We 
have performed numerical simulations of the disruption of a 
system like the Sgr dwarf orbiting a Galactic potential with 3 
components: a bulge, a disk and a dark logarithmic halo 
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where d — 12 kpc and Whaio — 13L5 km/s (Johnston et al. 
1996). The parameter q is allowed to vary from 0.8 to 1.25, 
that is, from an oblate to a prolate configuration^. The satel- 
lite galaxy is modeled by a set of 50000 self-gravitating parti- 
cles. We use a quadrupole expansion of the internal potential 
of the system, and choose a King profile for the pre-disruption 
dwarf. The orbital initial conditions are chosen to satisfy the 
constraints given by the present position and radial velocity of 
the main body of Sgr (Ibata et al. 1997). For each of the q 
values of the dark halo potentials, we select orbits which have 
similar (mean) pericenter and apocenter distances as well as 
comparable L^ (z-component of the angular momentum) satis- 
fying the above mentioned constraints. Our models of the dwarf 
are also slightly readjusted for each of the dark halo shapes, so 
as to produce the same remnant system by the present day. This 
implies that differences in the characteristics of the debris may 
only be attributed to a change in the flattening of the poten- 
tial. The other possible free variables, such as the model of the 
dwarf and the orbital parameters, are (essentially) the same by 
construction. 

We have studied five flattenings for the dark halo poten- 
tial q ~ [0.8, 0.9, 1.0, 1.11, 1.25], which correspond to average 
density axis ratios (within the orbit of Sgr) (qp) ^ [3/5, 4/5, 1, 
5/4, 5/3]. In this Letter we focus on the properties of the satel- 
lite debris after 10 Gyr of evolution. 

2.1. Young streams 

We study now the characteristics of young debris streams. 
We select those particles that have been released by the dwarf 
in the last three pericentric passages (less than 2 Gyr ago in our 
models). The panels in Figure 1 show the heliocentric radial 
velocity vs. longitude Aguir^ for this set of particles and for the 
different shaped dark halos. The large gray symbols in each 
of the panels correspond to the stars in Table 3 of M04 with 
distances is larger than 13 kpc. There appears to be very good 

^It is worth noting that the isodensity surfaces are more oblate (prolate) than 
the equipotential suifaces. 

'This variable corresponds to the angular distance along the stream. Its 
relation to (/, b) has been derived in Majewski et al. (2003) 
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Fig. 1 . — Distribution of (Asun, K-) for particles after 10 Gyr of evolution, for each one of the simulations with different degrees of flattening for the Galactic 
dark matter halo (black dots). Particles that have been released in the last 2 Gyr are shown in black, while the gray dots corresponds to those lost between 2 and 4 
Gyr ago. The large gray symbols denote the stars observed by Majewski et al. (2004) with distances larger than 1 3 kpc, which is the criterion used by these authors 
to minimize the contamination by foreground thick disk and halo stars. Solid symbols correspond to stars which have been considered by M04 as likely members of 
the streams, while open symbols denote those that have been statistically identified as outliers by these authors. The (random) velocity errors quoted in Majewski et 
al. 2004 are ~ 6 km/s. 



agreement between the models and the stars in the streams of 
Sgr observed by M04. 

To quantify how good the agreement is, we perform a x^ test. 
We bin the particles in the simulations (again only those lost in 
the last 3 pericentric passages) in elements of AAgun ~ 5° 
width, and measure the mean radial velocity and the disper- 
sion in each bin. We proceed in a similar way for the stars in 
M04 sample, selecting those that are more likely members of 
the stream as derived by M04 (these are the solid gray points in 
Figure 1). We then measure x^ as 
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where Vi is the mean radial velocity in the ith bin, and (Ji its 
dispersion. The subscripts o and m correspond to the observa- 
tional data and to the model points, respectively. The results of 
this test are shown in Figure 2 and in the last column of Table 1 . 

Figure 2 confirms that in all cases, the models provide a 
good representation of the data: The value of the reduced x^ 
(which we defined as x^/^bms) is close to unity. The most 
favored case (taken as that with the smallest x^) seems to be 
the most flattened halo q = 0.8. A bin to bin analysis of Fig- 
ure 2 shows that the bin located at Agun = 52.5°, which con- 
tains only one data point, produces the largest deviation. If this 
star is removed, the x^ is significantly reduced for all shapes. 
In this case, the most favored models are those with q = 0.8 
({c/a)p ^ 3/5) and q = 1.25 ({c/a)p ^ 5/3). These corre- 
spond to the most oblate and to the most prolate halos we have 
considered. 

We have also measured the x^ obtained by comparing mod- 
els to models for the same range of longitude Agun (see Ta- 
ble 1). We find that x^ is typically of the order of unity, and 
that the range in x^ values is the same as that obtained in the 
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Fig. 2. — The predicted (asterisks) and observed (gray diamonds) distribu- 
tion of mean heliocentric radial velocity measured in bins of AAsun = 5° 
width, for the particles (in black) of Fig. 1 and for the stars shown as filled gray 
symbols in the same figure. 



data-models comparison. Small values of x^ are found for the 
comparison of model q — 0.8 with q — 0.9 and with q = 1.25, 
as well as for the comparison of q = 1.00 with q = 1.11. 
This is most likely due to the fact that the sky projected orbital 
paths of the center of mass of Sgr overlap almost perfectly up 
to Asun ^ 100° for each of these two groups of models. 

The extent of the stellar stream observed by M04 can be com- 
pared to our simulations to derive its age, i.e. the time elapsed 
since the stars became unbound from the Sgr dwarf. Since mass 
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loss preferentially occurs during pericentric passages, the age 
of a (portion of a) stream may also be measured by the num- 
ber of pericentric passages that have passed since its formation. 
This quantity can easily be derived for the streams in our sim- 
ulations. The material found in the range Agun = [20°, 150°] 
has been lost in the last 3 passages: ^0.1 Gyr ago (domi- 
nates close to the bound core), ^ 0.85 Gyr ago and ~ 1.6 Gyr 
ago (dominant at Asun<;100°). Uncertainties in this estimate 
can be attributed to incomplete knowledge of the orbit of the 
dwarf or of its mass-loss rate. The former, however, is well 
constrained by the available data (its average radial period lies 
in the range 0.7 - 0.76 in our models). The mass-loss rates 
are, on the other hand, constrained by the velocity dispersion 
found in the streams. Higher mass-loss rates than those found 
in these simulations would lead to thicker/hotter streams and at 
a given location, the stream will be younger. On the other hand, 
lower mass-loss rates would lead to colder structures (since the 
particles lost would have a smaller range of energies) which 
would be older on average. Inspection of Fig. 1 suggests that 
the uncertainty in the quoted values cannot be larger than ±1 
pericentric passage, that is ±0.75 Gyr. 

2.2. Older streams 

How old do streams have to be before they show noticeable 
dissimilarities due to variation in the Galactic dark halo shape? 
The kinematics of particles lost between 2 and 4 Gyr ago (gray 
dots in Fig. 1) already exhibit important and measurable differ- 
ences. These are further highlighted in Figure 3, which plots 
the trend in the heliocentric radial velocity as a function of lon- 
gitude Asun for three different cases: q = 0.8 (light gray; dia- 
monds), q — 1.00 (dark gray; asterisks) and q = 1.25 (black; 
triangles) for particles lost in the last 4 Gyr. The error bars 
indicate the expected velocity dispersion around the mean he- 
liocentric radial velocity in a given longitude bin. The different 
models are clearly distinguishable at Agun ^ 300° in the trail- 
ing stream, and for Agun ~ 240° in the leading stream. 

In fact Law et al. (2003) show in their Figure la, the radial 
velocities of a set of M giants from the 2MASS catalog located 
in the leading stream at Agun = [200°, 360°]. We have dig- 
itized this plot (since the data are not yet available in tabular 
form) and show the corresponding data points as open squares 
in the top panel of Figure 4. The mean radial velocities mea- 
sured as function of Agun for the most oblate, the spherical and 
the most prolate models are overplotted for comparison. In or- 
der to reproduce the behavior of the radial velocity data around 
Asun ^ 240° in the leading stream, a prolate halo is required. 
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Table 1 
Reduced x^ for the different experiments. The first 4 

COLUMNS are obtained BY COMPARING THE DIFFERENT 

MODELS TO ONE ANOTHER. THE LAST COLUMN CORRESPONDS 

TO THE DATA TO MODEL xled- THE VALUE BETWEEN BRACKETS 

IS OBTAINED WHEN THE BIN AT Asun = 52.5° IS NOT 

CONSIDERED. 



A clear trend exists in the sense that the smaller the axis ratio of 
the dark halo, the more negative the mean galactocentric radial 
velocity at this location is, and hence, the less likely the model 
becomes. If the data in Law et al. (2003) are confirmed, this 
would be the first direct evidence supporting a prolate shape 
for the dark halo of our Galaxy. 

While the q — 1.25 model gives the best match to this new 
dataset, it seems clear that an even more elongated halo could 
yield a better fit. Another option could be to reduce the mass of 
the disk by a small amount while keeping this axis ratio fixed. 
We explore these two possibilities using test particle orbits. The 
results are presented in the bottom panel of Figure 4 for the lon- 
gitudes probed by the leading and trailing streams traced so far. 
The black curves correspond to our original set of parameters: 
Mdisk = 10"Mo, Uhaio = 131.5 km/s, and rf == 12 kpc, for 
the cases q ^ \ (dashed), q — 1.25 (solid) and q ~ \.b (dot- 
ted). The dark gray curves have been computed with a lower 
disk mass Mdisk = 5 x 10^°Mq, while keeping the vlsr and 
the circular velocity at 200 kpc fixed (this changes Uhaio ~ 137 
km/s and d = 3.5 kpc). Therefore we conclude that reducing 
the mass of the disk has a negligible effect on the orbital path. 
Moreover, the very prolate q — l.h halo seems to be too elon- 
gated and does not trace well the trends in the data shown in the 
top panel of Figure 4. 

3. SUMMARY 

There is good agreement between the models and the stars 
in the trailing stream of Sgr observed by M04. The compar- 
ison confirms that the different models can scarcely be dis- 
tinguished from one another when only the dynamics of the 
youngest streams (< 1.6 Gyr of age) are analyzed, since all 
of them reproduce to similar extent the radial velocity trends 
observed in the data. 

To make reliable estimates of the shape of the dark halo of 
the Galaxy slightly older streams are needed. The Sgr leading 
stream data published by Law et al. (2003) encompass such 
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Fig. 3. — The predicted mean heliocentric radial velocity measured in bins 
of AAsun = 5° width, for particles lost in the last 4 Gyr for models with 
q = 0.8 (light gray, diamonds), q = 1.00 (dark gray, asterisks) and q = 1.25 
(black, triangles). The solid curves indicate the actual trajectory followed by 
the center of mass in the simulations. Note that while the overall behavior of 
the latter is similar, the actual mean value at each longitude bin can be quite 
different. The error bars indicate the velocity dispersion around the mean in 
each bin. 



A. Helmi 





; q=0.80 




3 = 


1.00 


q^ 


1.25 \ 


200 


■"XSa 












^ 


100 


^ ^ 










i. 


C \ 





\ 








\ aw? 


¥ 


^ 


100 


[ 


^. V 


"SC Q„ 




^^^wW' 




\ 


200 


: (o) 


\ 










H 












--^ 




I 


















200 


■~^ 








'' /^ 


% 


X 1 


100 


r A 


V 








/ 


^1 





: 


\ 






,;'^ / 


/ 


' 1 


100 


- 


\ 


\^ 


^ 


^I^^^Va*^ 




1 


200 


: (b) 












1 



300 200 100 

Longitude A^^^ (degrees) 



Fig. 4. — Panel (a): The predicted mean galactocentric radial velocity vs 
Asun measured in bins of AAsun = 5° width, for particles lost in the last 4 
Gyr of evolution for models with q = 0.8 (light gray; dot-dashed), q = 1.00 
(dark gray; dashed) and q = 1.25 (black; solid). The asterisks correspond 
to the data published by M04, while the open squares are from Law et al. 
(2003). Panel (b): The galactocentric radial velocity vs Asun for test particle 
orbits integrated in different potentials: q = 1 (dashed), q = 1.25 (solid) and 
q = 1.5 (dotted). The black curves correspond to the original parameters, 
while the dark gray to halving the disk mass. 



older debris (^^ 2 — 4 Gyr of dynamical age). This dataset 
strongly suggests that the dark halo of the Milky Way is prolate 
and that the average density axis ratio {c/a)p within the orbit 
of Sgr is probably larger than 5/3 (q = 1.25) but smaller than 
12/5 (q = 1.5). 

This result has a range of important implications. It rules out 
models in which the dark matter is distributed as the baryons in 
our Galaxy, such as for example in the form of cold H2 clumps 
in the outer disk (e.g. Pfenninger & Combes, 1994). Moreover, 
it poses a big problem for MOND because requires the MOND 
force field - which depends on the (oblate) distribution of lu- 
minous matter in the Galaxy - to mimic the gravitational effect 
of a prolate halo. Finally, a prolate shape could go a long way 
in explaining the Holmberg effect (e.g. Sales & Garcia-Lambas 
2004): the excess of satellites observed along the minor axis of 
disk galaxies may just be reflecting the shape of the underlying 
mass distribution. 

I would like to thank Heather Morrison, Paul Harding, Ed Ol- 
szewski, Simon White and Mariano Mendez for very enjoyable 
discussions, and the anonymous referee for a very constructive 
report. 
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